Nodeless d-wave superconductivity in a quasi-one-dimensional 
organic superconductor under anion order 



On 
On 
On 

o 

Q 

00 



S3 
O 

o 

I 1 
5h ■ 

33 ; 



I 

o 
o 



> 
m 

(N 

On 
OS 
-t— > 



i 

O 

o 



Hiroshi Shimahara 

Department of Quantum Matter Science, ADSM, Hiroshima University, Higashi- Hiroshima 739-8526, Japan 

(Received December 1999) 

We propose a mechanism of nodeless d-wave superconductivity in a quasi-one-dimensional organic 
superconductor under anion order. It is shown that split of the Fermi surfaces due to the anion order 
eliminates the line nodes of a d-wave gap function on the Fermi surfaces, and a small finite energy 
gap appears in the quasi-particle excitations of superconductivity. We discuss that temperature 
dependences of thermal conductivity and NMR relaxation rate observed in a quasi-one-dimensional 
organic superconductor (TMTSF^CICU can be explained consistently by the small energy gap. 

PACS numbers: 74.70.Kn, 74.20. Mn, 74.20.-z 

Anisotropy of superconducting gap in organic super- 
conductors is an important subject, since it is closely re- 
lated to the pairing mechanism. In particular, a mecha- 
nism of an anisotropic superconductivity by pairing inter- 
actions mediated by antiferromagnetic (AF) fluctuations 
was proposed by Emery in the quasi-one-dimensional 
(Q1D) superconductors (TMTSF) 2 X §, and has been 
studied by many authors, because of the proximity to the 
spin density wave. For example, it was shown that the 
phase diagrams and line nodes are semi-quantitatively 
reproduced by this mechanism in a microscopic calcula- 
tion B. At the present, however, there is not any exper- 
imental evidence that clarifies the origin of the pairing 
interactions in those compounds. 

For this subject, an NMR experiment was made in 
a Q1D superconductor (TMTSF) 2 C10 4 by Takigawa et 
al. |J. In their data, the nuclear relaxation rate T^ 1 did 
not exhibit a coherence peak (Hebel-Slichter peak) just 
below the superconducting transition temperatue T c , and 
it was proportional to T 3 at lower temperatures. These 
results suggest that the superconducting gap function is 
anisotropic and has line nodes on the Fermi surface. By 
calculating the relaxation rate, Hasegawa et al. proposed 
two possibilities, singlet pairing with the gap function 
A(k) ~ cos(fcj,) + C (where C < 1) and triplet pairing 
with A(k) ~ sm(k y ) §). Here, we choose x-axis in the 
direction of the highly conductive chain, i.e., a-axis, and 
the y-axis in the direction of &-axis. 

However, recently, Belin et al. have observed a tem- 
perature dependence of thermal conductivity which indi- 
cates the nodeless gap |)|]. It may appear to contradict 
the NMR data at low temperatures, but as they pointed 
out, the NMR experiment was made only for T > 0.6 T c , 
which may not be sufficiently low to determine the power 
of T. On the other hand, the absence of the peak of T^ 1 
just below T c still suggests an anisotropic superconduc- 
tivity. Hence, they argued the possibility of the triplet 
superconductivity with a gap function A(k) ~ sin(fc x ), 
which is anisotropic, but is finite over the whole open 
Fermi surfaces. 

For this triplet pairing, the Hebel-Slichter peak is much 
smaller than that for the s-wave pairing since the coher- 



ence factor is canceled in the momentum integral which 
gives the relaxation rate. However, in a theoretical pre- 
diction , it does not seem to be small enough to explain 
the experimental data in which the peak is not seen at all. 
Further, although the NMR experiment was not made at 
low temperatures, a tendancy of T^ 1 oc T 3 was seen 
above 0.6 T c . Thus, if we believe the NMR data at the 
present, it is plausible that an anisotropic gap function 
which gives a quasi-particle density of states similar to 
that given by a gap function with line nodes developes 
at least for T > 0.6 T c . 

In addition to the inexplicability in the NMR data, 
the following two points are not clarified in explanations 
based on the triplet pairing: (1) The triplet pairing seems 
to contradict that the proximity to the antiferromag- 
netism may favor the <i-wave superconductivity; (2) If 
we assume the triplet pairing, a reentrant superconduct- 
ing phase must be observed at a higher field above the 
upper critical field as Lebed' and Dupuis et al. calcu- 
lated , but there is not any experimental evidence of 
such reentrant phase at the present g] . 

In this paper, we propose a mechanism of nodeless d- 
wave superconductivity which could explain the behav- 
iors mentioned above, which appear to be inconsistent 
with each other. First, we assume pairing interactions 
mediated by the AF fluctuations. In the absence of an- 
ion order, the line nodes are obtained near k y — ±7r/2, 
as in a result of a microscopic calculation When an- 
ion order occurs, the lattice periodicity changes in the 
6-direction, and the Brillouin zone becomes half in the 
fc a -direction. Since the degeneracies are removed by po- 
tentials due to the anions, which are staggered in the b- 
direction, the Fermi surfaces split at the new zone bound- 
ary, near which the line nodes are situated. Hence, the 
line nodes are eliminated for appropriate parameters. 

Recently, Yoshino et al. estimated the staggered po- 
tentials due to the anion order E g w 0.083 1 from their 
experimental data of magnetoresistance where t de- 
notes the hopping energy along the chain. We will show 
below that the value of E g is, as an order of magnitude, 
large enough to eliminate the nodes of the cJ-wave super- 
conducting gap near k y — ±tt/2. 
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We start with a two-dimensional tight-binding Hamil- 
tonian 
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H = 



E 
E 

1,3 



(1) 



with effective interactions Vij. We define two sub-chains 
A and B, and site energies e,; = £ s and — E g for i E A 
and i e B, respectively. The hopping energies i»j are 
assumed to be zero unless the sites i and j are nearest 
neighbor. For the nearest neighbor sites, ty = t when i 
and j are on the same chain, while Uj = t' when i and j 
are on different sub-chains. We use units with t = 1. 

The bilinear terms of the Hamiltonian are diagonalized 
by the transformations 



E 



» 



(2) 



with a = a, 6 and s 



where we define a, 



for i e A and bj t „ — c JiCr for j 6 £?. We obtain 
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1/2 



(3) 



U a - = U b+ = V{k y ) = f* 

and the dispersion relations of the (s) electron bands 

e«(k) = -2tcos(fc x ) - sE y (k y ), (4) 



where = y £ y 2 + E g 2 and e y = —2t'cos(k y ). These 
forms of the dispersion relations were used by Yoshino et 
al. to explain their experimental results |^]. The Fermi 
surfaces are modified by the staggered potential as shown 
in Fig.|l[ In particular, the Fermi surfaces split at the 
boundary of the new Brillouin zone (k y — ±7r/2). 



We write the Fourier transformation of Vij as 
V(c l ) = V (q x ) + Vi(< l ) 



(5) 



with intrachain and interchain interactions Vq and V\ 
respectively. Since the spin fluctuations are enhanced 
by the Fermi surface nesting with nesting vector q TO « 
(±7r, ±7r/2), which is equivalent to q m sa (0, ±7r/2) in 
the presence of the anion order, the spin susceptibility 
x(q, cj) has peaks at q = q m . The peak width is very 
small in the p^-direction, while it is broad in the p y - 
direction (|]. Therefore, the effective interactions due to 
the exchange of the spin fluctuations are of long range in 
the a-direction, but of short range in the ^-direction jio) . 
Thus, we take interactions up to those between nearest 
neighbor chains 



Vi(q) = -vi(q x )cos(q y ). 



(0) 



The functions Vo(q x ) and Vi(q x ) have very sharp peaks 
at q x = ±7r/2, reflecting the behavior of x(q, 
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FIG. 1. Fermi surfaces of the (+) and (— ) electron bands 
for 1/4 filling, when t' = 0.1. Solid and broken lines show 
those for E g = 0.083 and E g — 0, respectively. 

Here, we make a simplification that the peak widths 
of Vo(q x ) and vi(q x ) are similar, which is qualitatively 
justified when t' /t ^ 1. Further, in order to express the 
sharp peak, we employ Lorentzian form 



E 

^=±tt/2 
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(7) 



and put V (q x ) = gof(q x ) and V\(q x ) = gif(q x )- 

Since the peaks are situated at the nesting vectors 
which connect the (+) and (— ) electron bands, the peak 
width corresponds to the width of the area near the Fermi 
surfaces in which the pairing interactions are efficient. 
Therefore, we can assume that important contributions 
to the superconductivity mainly come from the pairing of 
electrons near the Fermi surfaces as in the starndard weak 
coupling theory. Hence, we can put (c^t cS^ ) = 
when Eg and t' are not too small, since k and — k can 
not be near the Fermi surfaces of (=f) and (±) electron 
bands, respectively, at the same time. This approxima- 
tion is appropriate qualitatively when E g 3> k^T, and 
the experimental result E g = 0.083 1 in (TMTSF) 2 C10 4 
satisfies this condition. 

Hence, we obtain a gap equation 

A( s )(k) = -■^5^V("'J(k,k / )W('''(k / )A('')(k / ) ) (8) 

where we define 
jy-(s')(k') 
£( s ')(k') 



k's' 



2£( s ')(k') 



, ,£ (s,) (k'K 
— tanh( — ) 



2k B T 



(9) 



V(e( s ')(k')) 2 + (A( s ')(k')) 5 
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The effective pairing interactions are obtained as 



1 



\/(+ + )(k,k') = V { 
= i[vb(k-k' 

F(+-)(k,k') = V { 
= ±[vb(k-k')(l 



>(k,k') 

I E 9~ 

E v E 'v 

)(k,k') 

E„ 2 



+ ^l(k-k') J 



V^(k-k') 1 



where E' 



(10) 



. y - Eyik'y) and e y ' = e y (k' y ). 
Since the contributions from the electrons near the 
Fermi surfaces are dominant as we discussed above, we 
replace the effective pairing interactions with those on 
the Fermi surfaces. Further, for the sharp peak of the 
pairing interactions, we introduce an effective cutoff en- 
ergy lo c . Thus, we rewrite the right hand side of eq.(8) 
as 



^V*>(o,fc;) / dew^\K,k' y ) 

-tt/2 K J-u> c 

xK^(k y ,k')A^'\k'), 



(11) 



where we define 



K^(k y ,k' y ) ee V^'H(k x ,ky),(k&(K)> k 'y)) 

+ V^'H(k X ,ky),(-k { F i ) (k>),k>)) 



A^(ky) EE A^(±k^(ky),ky). 



(12) 

Here, the function k^(k y ) > gives the value of k x on 
the Fermi surface of the (s) electron band at k y . The 
function p^ s \0, k y ) expresses the density of states at e — 
and k y of the (s) electron band of k x > 0, which is 

calculated as p^(0,k y ) = l/|47rt sin(fc^ ) (fc y ))|. 

Therefore, in the limit of T — > T c , we obtain the lin- 
earized gap equations in a form of eigen value equations 



AAM(Av) 



*/ 2 dk'„ 



XK^'Hky,k' y )M S '\k' y ), 



(13) 



where A = 1/ log(2e 7 o; c /7rfcBT) with the Euler constant 
7 = 0.57721 • • •. The eigen solution which belongs to 
the maximum positive eigen value gives the momentum 
dependence of the gap function near T = T c . 

We solve eq.(13) numerically for E g — 0, 0.083, and 
0.2, when t' = 0.1. The chemical potential p is adjusted 
so that the filling of holes n = 1/4. Since it was found 
in a microscopic calculation j|] that the range of the in- 
teraction are over 30 sites along the chain, and that in- 
terchain interactions are much weaker than intrachain 
interactions, we use the parameters q^jn = 0.01 and 
9o/9i = 5 as an example. The qualitative result that 
we will show below is not very sensitive to the values of 
these parameters. 



Figure g shows the numerical result of A' 
the normalization constant A is defined by 



A 



^ 2 dk 



^\A(ky)\ 



-tt/2 



1/2 



\k y ), where 



(14) 



which must increase as temperature decreases. When 
Eg = 0, the result in the absence of anion order is re- 
produced. We have usual ci-wave superconductivity with 
line nodes at k y « ±7r/2 (dotted lines). The gap func- 
tions A( + ) and A^ - ) have opposite signs equivalently to 
that A(\k y \ < n/2) > and A(vr/2 < \k y \ < tt) < 0, 
which corresponds to the pairing with A(k) ~ cos(k y ) in 
a more simplified model. When E g 7^ 0, we find that the 
line nodes vanish but the gap functions A' + J and A^' 
still have opposite signs. The magnitudes lA^I take 
minimum values at the zone boundary, and it is found 
that the minimum values increase with E a . 
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FIG. 2. The gap functions on the Fermi surfaces. The 
solid, dotted broken, and dotted lines show the solutions for 
E g = 0.2, 0.083, and 0, respectively. 

Figures || and [| show the densities of states near the 
Fermi level for E g = and 0.2, respectively. It is found 
that the energy gap opens when E g ^ 0, but except the 
small energy gap near E = the behavior of the density 
of states is very similar to that for E g = 0. 

The eigen value A can be regarded as an effective di- 
mensionless coupling constant. The value of Xt/go is es- 
timated to be 29.9, 25.1, and 14.8, for E g = 0, 0.083 1, 
and 0.2 1, respectively. One could make an estimation of 
T c if the values of go ft and to c are given, but in practice 
since an effect of self-energy renormalization, especially 
that from the pseudo gap, is very strong [0, we need a 
more careful treatment for a quantitative discussion. 
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FIG. 3. Density of states near the Fermi level, when 
Eg — 0. Solid and broken lines show the results for the (+) 
and ( — ) electron bands, respectively. 
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FIG. 4. Density of states near the Fermi level, when 
Eg — 0.2. Solid and broken lines show the results for the 
(+) and (— ) electron bands, respectively. 

In conclusion, we propose a mechanism of a nodeless 
d-wave superconductivity under the anion order. The 
momentum dependence of the gap function is calculated 
in a model with pairing interactions mediated by the AF 
fluctuations. It is found that the line nodes are located 
near the boundary of the half Brillouin zone in the ab- 
sence of the anion order, and because of their locations, 
the line nodes vanish in the presence of the anion order. 

Thus, the energy gap of the quasi-particle excitations 
of superconductivity becomes finite as shown in Fig||. 
Such a small energy gap would not be observed very well 
in experiments at high temperatures near T c . Since the 
behavior of the density of states is very similar to that for 
Eg = except the small energy gap, and the gap func- 
tion changes its sign depending on s = ±, temperature 



dependences near T c of observed quantities are expected 
to be similar to those for the d-wave gap function with 
line nodes. In particular, it is expected that a Hebel- 
Slichter peak does not occur just below T c and a ten- 
dancy of T£"~ oc T 3 is observed near and below T c , in 
NMR relaxation rate. 

On the other hand, we could expect behaviors of a full 
gap superconductivity, for low temperatures. As temper- 
ature decreases, the superconducting gap developes, and 
finally the temperature becomes smaller than the finite 
energy gap. For such low temperatures, the electronic 
thermal conductivity should become vanishingly small. 

These behaviors may explain the recent experimen- 
tal results of the thermal conductivity and the NMR 
relaxation rate in (TMTSF) 2 C104 consistently. How- 
ever, quantitative confirmation by fitting the experimen- 
tal data seems to be difficult, because detailed informa- 
tions of the pairing interactions, such as those in mo- 
mentum, frequency, and temperature dependences, are 
needed for this purpose. In the future, observation of the 
Knight shift may clarify whether the pairing is singlet or 
triplet. 

This work was supported by a grant for CREST from 
JST. 
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